Analyzing flow anisotropies with excursion sets in relativistic heavy-ion colhsions 
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We show that flow anisotropies in relativistic heavy-ion collisions can be analyzed using a certain 
technique of shape analysis of excursion sets recently proposed by us for CMBR fluctuations to 
investigate anisotropic expansion history of the universe. The technique analyzes shapes (sizes) of 
patches above (below) certain threshold value for transverse energy/particle number (the excur- 
sion sets) as a function of the azimuthal angle and rapidity. Modeling flow by imparting extra 
anisotropic momentum to the momentum distribution of particles from HIJING, we compare the 
resulting distributions for excursion sets at two different azimuthal angles. Angles with maximum 
«) ' difference in the two distributions identify the event plane, and the magnitude of difference in the 

, two distributions relates to the magnitude of momentum anisotropy, i.e. elliptic flow. 
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Analysis of anisotropies of flow in relativistic heavy-ion collision experiments (RHICE) has provided valuable in- 
formation regarding the properties of matter produced during early transient stages in these collisions 1]. Early 
investigations focused on the elliptic flow relating to the 2nd Fourier coefficient of momentum anisotropy as a function 
of azimuthal angle, as well as couple of other flow coefficients. Drawing correspondence of the underlying physics of 
fluctuations between RHICE and the Cosmic Microwave Background Radiation (CMBR) anisotropies, it was pointed 
out by us in ref. 0] that a wealth of information is contained in plots of root-mean-square values of flow coefficients 
^ ' f™* for a large range of values of n, from 1 to 30-40. It was also proposed in [3] that the presence of superhorizon 
^ fluctuations can lead to the suppression of superhorizon modes, as well as the existence of acoustic oscillations in plots 
of v™" in RHICE, similar to the celebrated acoustic peaks of CMBR arising from inflationary density fluctuations. 
I This suppression of superhorizon modes in RHICE may have been seen in the analysis reported by Sorenson for RHIC 
> ' data [3|]. 

t — ' This very useful interplay of physics from CMBR to RHICE was further explored in ref. [i!] following the effects of 
magnetic fleld on the power spectrum of CMBR anisotropies Utilizing that insight, it was proposed in Q that 
' presence of strong magnetic field in the initial stages of a non-central collision in RHICE can affect flow anisotropies, 
[ , for example, it can lead to enhancement of flow coefficients. In particular, it could lead to enhancement in the elliptic 
flow coefficient V2 by almost 30 % (which may allow for a higher value of viscosity) . A different analysis in ref. [6] in 
the strong field limit further explored these effects and gave similar results for the enhancement of V2 ■ 

In this paper we explore another common feature between the analysis of CMBR and the flow analysis in RHICE. 
. . \ This is based on a technique recently proposed by us for investigating any anisotropic expansion stage in the history 
of the universe by analyzing the shapes of patches of CMBR fluctuations (the excursion sets) at the surface of last 
scattering Q. The main idea underlying this analysis Q was the following. If the density perturbations generated 
initially are statistically isotropic, then they will get deformed in a specific direction if the universe ever went through 
' an anisotropic expansion. This (statistically averaged) deformation can be detected by a shape analysis of fluctuation 
patches of CMBR. An anisotropic expansion of the universe has very strong similarities with the anisotropic flow 
expected in RHICE. The situation of anisotropic expansion of the universe is directly realized in the form of elliptic 
flow for non-central events. Flow anisotropies are generally present in every event in RHICE and all values of v™" 
are non-zero in general due to the presence of fluctuations of different scales even in the central collisions 3. This 
suggests that the shape analysis of excursion sets used in ref. [7|] should provide useful information for flow anisotropies 
in RHICE as well, especially for the elliptic flow. We will see below that this is indeed true. We consider fluctuations of 
transverse energy or particle number as a function of the azimuthal angle and rapidity. This is like CMBR fluctuations 
at the last scattering surface. Then by considering fluctuations above/below a certain threshold value, we generate 
the excursion sets of fluctuations, we call them simply as fluctuation patches. To see the possibility of detecting 
anisotropic flow, we model flow by imparting extra anisotropic momentum to the momentum distribution of particles 
obtained from HIJING ^8J. We then analyze the shape/size distribution of the fluctuation patches at two different 
azimuthal angle. We will see that azimuthal angles with maximum difference in the two distributions identify the 
event plane, and the magnitude of difference relates to the magnitude of momentum anisotropy, i.e. elliptic flow. 
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FIG. 1: 

Plot of excursion sets (fluctuation patches) in the azimuthal angle (j> (in radians) and the rapidity y space. Horizontal solid 

lines show the slicing used for widths of filled patches. 

For analyzing anisotropic expansion history of the universe in [7| , we had examined the size distribution of patches 
on the surface of last scattering (which is a two sphere S'^) along 9 and (j> directions. This will be different for RHICE. 
Since our interest is in probing flow anisotropy, say the elliptic flow in a non-central event, we would like to compare 
the shape/size distribution of fluctuations along the X axis with that along the Y axis, where the X and the Y axis 
correspond to the event plane. As we mentioned, we model anisotropic flow in HIJING by multiplying a factor fp^ 
to the momentum of every final state particle. For central events fp^ is taken to be 1 while for non-central events 
fp^ is taken to vary from 1.1 to 1.4 representing the momentum anisotropy resulting from elliptic flow. Precise value 
of momentum anisotropy represented by fp^ is not an issue here (though the values used here are of right order of 
magnitude, say, fp^ = 1.2 implying about 20 % anisotropy). Our intention is only to show the systematic pattern in 
which the difference between the size distributions of excursion sets at different azimuthal angle probe the event plane 
and strength of elliptic flow. A quantitative correspondence can only be established by incorporating our analysis in 
hydrodynamical simulations and we hope to report such an analysis in a future work. 

With the (anisotropic) momentum distribution of particles from HIJING, we collect particle numbers as well as 
transverse energy in different (p and y bins where is the azimuthal angle in the event plane and y is the rapidity. 
We use 100 bins for full 27r range of (j> and 100 bins for the full range of y. We have also considered smaller range of 
rapidity y and the results are similar. Full range of y gives better statistics. To exclude very large values of y we use 
a lower cutoff for Pt — 150 MeV. For each event thus we get the particle number/transverse energy as a function of 
(j) and y. By subtracting the average values of these quantities we get particle number fluctuations/transverse energy 
fluctuations as a function of (p E^nd rj. We then consider a threshold value for these fluctuations and consider only 
those (j), y bins where particle number or transverse energy fluctuations is above (or below) this threshold value. This 
provides us with the excursion sets for the relevant fluctuations. Fig.l shows the excursion sets for transverse energy 
fluctuations obtained from HIJING for a single Au-Au collision event at ^/snn = 2 TeV. (Here, we will be presenting 
results for transverse energy fluctuations, results for number fluctuations are very similar to these, hence we do not 
discuss them separately.) We consider here non-central events with the impact parameter lying between 6 = 3 — 4 
fm. Here we have used fp^ —1.2 (for modiflcation of the momenta of particles from HIJING) and fluctuations above 
0.2 X maximum fluctuation have been included to generate the excursion sets. We use here large collision energy as 
relevant for LHC. The size distributions for small number of particles have large statistical errors. For the case of 
RHIC energies one will need to combine much larger number of events to control statistical errors. 

What we want to know is the average width of the fluctuations along direction near cj) — and compare it with 
the similar average width near ip — tt/2. This will probe the elliptic flow. (Clearly, the same result will be obtained 
by comparing distributions at = tt and (f) ~ Sir/ 2.) We determine the distribution of widths along the (j) direction 
in a 45° wide arc centered at </) = and compare it with the distribution of widths (again, in a 45° wide arc) centered 
at (/) = 7r/2. To calculate these distributions we proceed as follows. We divide the entire region shown in Fig.l (with 
y varying within full range, and (j) ranging from to 360°) into thin slices (with a rapidity width of 0.18) along the 
direction, as shown by horizontal lines in Fig.l. Using these slices, we determine the (j) widths of various fllled patches. 
We then plot the frequency distributions (histograms) of the widths of the intersections of all the patches with these 
slices at = and at (j) — 'k/2 (within 45° arcs). For the isotropic expansion case, without any flow anisotropy, we 
expect these two histograms to almost overlap. For anisotropic expansion case we expect that fluctuation patches 
will be stretched more at ip = tt/2 than at cj) — 0, as system expands faster along ip = (and <p ~ tt) directions for 
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FIG. 2: 

(a) Plots of the histograms of the widths (in radians) of filled patches for fp^ — 1. Smooth curves join the points which are 
marked with corresponding (\/iV) error bars. Distribution for widths at (j) — 0(0 = 7r/2) are shown by solid (dashed) plots, 
(b) Plots for the case with momentum modification, with fp^ — 1.2. The inset in (b) shows the plots for larger widths. 



a non-central collision. Any relative difference, between the two histograms will, therefore, imply the presence of an 
anisotropy of expansion (apart from the possibility of any initial anisotropy, we will discuss it below) . 

Fig. 2a shows the frequency distributions (histograms) of the widths of the intersections of the slices with the 
filled patches (excursion sets) at = (solid curve) and (f> = n/2 (dashed curve) for the case with fp^ = 1 (i.e. no 
momentum anisotropy. The horizontal axis corresponds to the widths of the slices (in radians) along 0, with histogram 
bin having width of 0.05 radians. The vertical axis gives the frequency N of the occurrence of the respective widths 
in all the slicings of excursion sets (such as in Fig.l, but now for fp^ = 1). The error bars denote the statistical 
uncertainty of VTV for the frequency N in each bin. As we are comparing distributions for two different data sets 
(one at = 0, the other at cf) = 7r/2), we normalize the distributions, as well as the errors, with the total number of 
particles included for each data set. We combine histograms of 500 events in HIJING with the X axis of each event 
(with same event parameters such as the impact parameter, collision energy etc.) representing the minor axis of the 
ellipse representing the overlap of the two colliding nuclei. We can see from Fig. 2a that the two histograms, are almost 
overlapping. This is interesting as we are considering here non-central collisions (with 6 = 3 — 4 fm). This shows 
that the momentum distributions and all fluctuations obtained from HIJING (without any momentum modification, 
fp^ = 1) are isotropic. Any anisotropy will result from effects of particle interactions, primarily from collective 
flow effects, and we will represent it by using values of fp^ different from 1. We mention here that the peak near 0.1 
radians in Fig. 2 will represent fluctuations of size order 1 fm at the distance of about 10 fm representing the size of 
initial nuclei. As discussed in ref.Q, initial fluctuations are expected because of the nucleon size of order 1 fm (as 
well as quantum fluctuations in nucleon coordinates.) Such a peak thus may have useful information about the scale 
of fluctuations in the early stages of system evolution. 

We now consider an anisotropic momentum distribution by taking fp^ — 1.2 to modify particle momentum distri- 
bution from HIJING. (Note that this is just to model the effects of flow.) Fig. 2b shows the plots of distributions of 
widths for this case. We see that the distribution at = 7r/2 is smaller than the distribution at = (beyond the 
error bars) for small widths while the situation is reversed for larger widths (as shown by the inset in Fig. 2b). This is 
expected as the stretching of fluctuation patches will be more significant at <j) — tt/2 due to elliptic flow. Stretching 
makes fluctuation patches larger, hence the distribution for large widths is higher for (j> = 7r/2 than for (f> = 0. This 
is suitably compensated at smaller angles where </) = 7r/2 distribution falls below the (f> = distribution. We flnd 
that the ratio of the heights of the two peaks is directly related to the value of fp^ representing the effects of elliptic 
flow. For fp^ = 1.1, 1.2, 1.3, we find this ratio to be 1.3, 1.6, 2.0. Thus, the magnitude of elliptic fiow can be probed 
directly with the comparison of the distributions of widths of fluctuation patches at cj) — and at (p — Tr/2. Note that 
our modeling of flow in terms of a factor like fp^ is ad hoc. Thus we do not try to derive any functional relationship 
between this factor and the ratios of the heights of the two peaks above (or try to give error bars for these etc.). 
Our purpose here is to show specific patterns of these distributions depending on momentum anisotropy. We also 
mention here that in actual hydrodynamic flow, the effect of flow on the momentum of particles will depend on the Y 
coordinate of the particle inside the plasma (in the event plane), being less important for larger Y values, as well as on 
the distance from the center. To more accurately represent the Y dependence of fiow anisotropy, we have compared 
the (p = distribution using fp,^ = 1.3 with the (j> = Tr/2 distribution with /p^ = 1.1. The difference between the 
two distributions is found to be closer to the difference obtained for the constant value of fp^ — 1.2 for both the 
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FIG. 3: 

Plots of histograms (for 500 events each) at </!) = 0o (solid curve) and at (^o + 'r/2 (dashed curve) for 0o = (a), 7r/16 (b), n/8 

(c), and Tv/4 (d). 

distributions. Since we are only interested in comparative changes in the two distributions, for simplicity, we will 
continue to use a constant value of fp^ for determining the difference between the (j) = and (p — tt /2 distributions. 

Of course this requires identification of the event plane so that statistics of a large number of events can be added 
up for controlled statistical errors. We now show that our technique can provide a new way to identify the event 
plane. In Fig. 3 , we show a series of plots of widths distributions at = 0o (solid plot) and at = 0o + 7r/2 (dashed 
curve) for the angle 0o = 0, 7r/16, 7r/8, 7r/4. We note that the ratio of the two peaks (solid curve and dashed curve) 
is largest for 00 = which represents the event plane. Thus, one needs to simply plot the two width distributions 
at angles 7r/2 apart and keep changing the initial angle until the two distributions have maximum difference. This 
will identify the event plane. (Though we note that maximum overlap of the two distributions seems to happen 
at (/fQ = 7r/8 and not at 7r/4 as one would have expected. This could be due to the large arc size of 45° used for 
calculating the individual distributions.) We have also compared the two distributions by varying the separation of 
at which they are evaluated. When the separation is larger than about 50-60° the difference becomes almost similar 
up to 90° separation. With the physics of ellptic flow in view, we choose 90° separated distributions for comparison. 
Fig. 3 is shown with collection of data from 500 events each. One will need to do this on event- by- event basis to be 
able to identify the event plane in each case. Fig. 4 shows the same plots as in Fig. 3, but now for single events. We see 
that statistical errors are larger. For larger collision energies, statistical errors will be under better control. Also, one 
may arrive at the event plane by an iterative process. First one can roughly identify event planes for each event (as 
in Fig. 4) . Then a large number of events can be combined for each choice of 0o to maximize the difference between 
the two plots. This can be done by systematically changing the event plane for each event, until the all-event-sum 
histograms (as in Fig. 3) lead to maximum difference between the two histograms. 

In this paper we have shown that a shape/size analysis of excursion sets can lead to a completely new way of 
identifying the event plane in RHICE and to quantify the magnitude of flow anisotropies. This method is novel, and 
completely independent of other conventional methods. This technique can, thus, be used to cross-check the results 
obtained by other conventional techniques. Further, the entire distribution of sizes of excursion sets, as shown in 
Fig. 2, will contain valuable information about nature of fluctuations and that of anisotropic flow. Also the width 
distribution along the rapidity should provide information about the longitudinal scaling of fluctuations. These issues 
can be probed using hydrodynamic simulations. Work is underway along these line and we hope to present it in future. 
We mention here that anisotropic expansion of the universe was discussed in ref. ^7] using Fourier transform technique 
as well as direct shape analysis of excursion sets. In this paper we have used the direct shape analysis method for 
investigating flow anisotropy. It will be useful to also use Fourier transform technique for this purpose. Note that 
this technique will detect any effects present in the system which lead to anisotropy of momenta. For example the 
presence of initial strong magnetic field (as well as induced electric field) will lead to anisotropies [1, Q and should 
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FIG. 4: 

Series of plots as in Fig. 3, for single event case each. 



leave imprints on these distributions of widths. Even the presence of any initial anisotropies of fluctuations, such as for 
deformed nuclei, could be detected in this method. Important issue then will be to distinguish any initial anisotropy 
(as from a deformed nucleus) from anisotropy resulting from anisotropic flow. In this respect this technique of direct 
shape analysis seems particularly effective, in comparison to, e.g. the Fourier transform technique, as discussed in 
detail in ref.Ja . We hope to address these issues in a future work. We point out that other methods such as Minkowski 
fucntionals [s], analysis of genus statistics of excursion sets etc., have also been very effectively used for analyzing 
morphology of fluctuations in CMBR as well as in large scale structure in the universe, and it will be fruitful to apply 
these techniques also for analyzing the nature of fluctuations as well as flow anisotropies in RHICE. 

We are very grateful to Abhishek Atreya, Partha Bagchi, and Anjishnu Sarkar for useful discussions and comments. 
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